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ABSTRACT 

Context. Ultraluminous infrared galaxies are among the brightest objects in the local Universe. They are powered by strong star 
formation and/or an AGN. They are also likely to be the progenitors of elliptical galaxies. The study of the structure and kinematics 
of samples of local ULIRGs is necessary to understand the physical processes that these galaxies undergo, and their implications for 
our understanding of similar types of galaxies at high redshift. 

Aims. The goal of the project is to analyze the structure, dust distribution, ionization state, and kinematics of a representative sample of 
22 ULIRGs. The galaxies in the sample undergo different merger phases (they are evenly divided between pre- and post-coalescence 
systems) and ionization stages (27% Hn, 32% LINER, 18% Seyfert, and 23% mixed classifications) over a wide infrared luminosity 
range (1 1.8<L /R /L <12.6), with some galaxies of low-luminosity. The main aims of this paper are to present the sample and discuss 
the structure of the stellar and ionized gas components. 

Methods. Our study relies on the use of integral field optical spectroscopy data obtained with the INTEGRAL instrument at the 
William Herschel Telescope. 

Results. The structure of the ionized gas as traced by different emission lines has been studied and compared with that of the stellar 
continuum. We find structural variations between the gaseous and the stellar components, with offsets in the emission peaks positions 
of up to about 8 kpc. Young star formation (as traced by the Ha emission) is present in all regions of the galaxies. However, for 
64% of ULIRGs in an early interaction phase, the young star formation peak does not coincide with the stellar maxima. In contrast, 
galaxies undergoing advanced mergers have a Ha peak that is located in their nuclear regions. In three of the studied ULIRGs, hard 
ionizing photons traced by the [O m]/i5007 line excite extra-nuclear nebulae out to distances of about 7 kpc. These regions do not 
show bright stellar emission, but are rather dominated by nebular emission. These galaxies have nuclei classified as Seyfert in the 
literature. Approximately 40% of the pre-coalescence ULIRGs exhibit shifts between the peaks of their red continuum and that local 
to the [O i]/i6300 line. However, some of these peaks are associated with the secondary stellar nucleus of the system. In contrast, the 
emission in post-coalescence ULIRGs is concentrated towards the nuclei. These results imply that evolution caused by a merger is 
ocurring in the ionized gas structure of ULIRGs. 

Key words, infrared: galaxies - galaxies:interactions - active:galaxies - technique:spectroscopy 



1. Introduction 

Ultraluminous Infrared Galaxies (ULIRGs, 10 I2 L G < L bo i~ 
L/«[8 - 1000/im] < 10 13 L o ) are a key galaxy class character- 
ized by the emission of the bulk of their energy in the infrared. 
They are as luminous as QSOs and, in the local Universe, twice 
as numerous (see reviews by Sanders & Mirabel 1996; Lonsdale 
et al. 2006). The existence of galaxies with a strong infrared ex- 
cess was reported by Low & Kleinmann (1968) and Rieke & 
Low (1972), but the study of ULIRGs as a galaxy class started 

Send offprint requests to: M. Garcfa-Marin (maca@phl.uni-koeln.de) 
'Based on observations with the William Herschel Telescope oper- 
ated on the island of La Palma by the ING in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astroffsica de 
Canarias. Based also on observations with the NASA-ESA Hubble 
Space Telescope, obtained at the Space Telescope and Science Institute, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc. under NASA contract number NAS5-26555. 



with their detection in large numbers by the IRAS satellite (e.g., 
Soifer et al. 1984; Sanders et al. 1988). Subsequently, optical 
long-slit spectroscopy was used to investigate their main en- 
ergy source (e.g., Veilleux et al. 1995). It has been generally ac- 
cepted that ULIRGs are powered mainly by intense star forma- 
tion, whereas the presence and relative importance of an AGN 
was a matter of debate (Genzel et al. 1998; Veilleux, Kim & 
Sanders 1999; Risaliti et al. 2006; Farrah et al. 2007). However, 
IR studies suggest that the contribution of AGN to the bolomet- 
ric luminosity of the system may be relevant in 15-20% of cases 
(see Nardini et al. 2008 and Risaliti et al. 2006). 

ULIRGs are characterized by interaction/merger processes 
that strongly drive their dynamics. The merger triggers the star- 
burst activity (i.e., the ultra-luminous phase), is responsible for 
the complex tidally-dominated morphologies (e.g., Borne et 
al. 2000; Bushouse et al. 2002; Farrah et al. 2001), dominates the 
ionized gas kinematics (e.g., Colina, Arribas & Monreal 2005), 
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and possibly transforms spiral galaxies into intermediate mass 
ellipticals (e.g., Genzel et al. 2001; Tacconi et al.2002; Dasyra 
et al. 2006). 

In a cosmological context, the importance of ULIRGs is sig- 
nificant because they appear to be the low-z analogs to at least 
part of the population of star-forming luminous galaxies giv- 
ing rise to the far-IR background. These are the sub-mm galax- 
ies (Smail, Ivison & Blain 1997) and the ULIRGs detected by 
Spitzer (Perez-Gonzalez et al. 2005). It has been shown that, in 
contrast to the situation in the local Universe, the importance 
of the contribution of ULIRGs and the less luminous LIRGs 
(Luminous Infrared Galaxies, 10 n L G < L^ /~ L/«[8 - 1000/vra] 
< 1O 12 L ) to the IR extragalactic light increases with redshift 
(Lagache, Puget & Dole 2005; Le Floc'h et al. 2004; Lonsdale 
et al. 2004; Yan et al. 2004; Caputi et al. 2006a, 2006b). They 
are the main contributors to the comoving star density of the 
Universe at z>1.0 (Elbaz et al.2002; Le Floc'h et al.2005; 
Perez-Gonzalez et al. 2005; Caputi et al. 2007). Generally, these 
high redshift studies do not focus on the detailed characteriza- 
tion of individual objects, but rather on their general integrated 
properties. To be able to study spatially resolved information, it 
is however important to understand their internal structure, dust 
distribution, ionization structure, and kinematics. This informa- 
tion is limited at high redshifts, and studies have used integral 
field spectrographs to analyze these high redshift objects in de- 
tail (e.g., Genzel et al 2006; Forster-Schreiber et al. 2006, 2009). 

In the local Universe, where the spatial resolution allows 
us to study the internal galaxy structure in greater detail, com- 
prehensive studies of local ULIRGs have been completed using 
imaging and long-slit spectroscopy in the optical and infrared 
(e.g., Veilleux et al. 1995; Genzel et al. 1998; Kim et al. 1998; 
Scoville et al.2000; Farrah et al.2001; Bushouse et al.2002; 
Imanishi et al. 2007). In some cases, the combination of these 
images with optical two-dimensional spectral information such 
as that provided by integral field spectroscopy (IFS) has been 
used to carry out comprehensive studies of individual or a few 
ULIRGs (e.g., Arribas, Colina & Borne 1999; Monreal-Ibero 
2004; Colina, Arribas & Monreal-Ibero 2005). 

To obtain statistical knowledge of the nature of ULIRGs, we 
have started a program to study, on the kpc scale, the dust struc- 
ture, the two-dimensional ionization statem, and the ionized gas 
kinematics of a representative sample of local ULIRGs, based 
on optical IFS data obtained with the instrument INTEGRAL. A 
similar study but for a representative sample of LIRGs is being 
conducted using VIMOS at the VLT (Arribas et al. 2008) and 
PMAS on the Calar Alto 3.5 m (Alonso-Herrero et al.2009). 
All these studies are part of the same project, which investi- 
gates the two-dimensional extinction, ionization, and kinematic 
kpc-scale structure of a representative sample of low-z LIRGs 
and ULIRGs, along with its implications for their high-redshift 
analogs. For instance, the IFS data presented in this paper sam- 
ple regions similar to those that will be studied by the instru- 
ments MIRI (Mid IR Instrument) and NIRSpec (Near Infrared 
Spectrograph), due to be launched onboard of the James Webb 
Space Telescope, at redshifts above 1 . 

In this paper, we combine data from previously published 
works (see Colina, Arribas & Monreal 2005; Garcfa-Mann et 
al. 2006 and references therein) and new data of 13 unpublished 
galaxies. The paper is organized as follows. We describe the 
sample selection in Sect. 2; the observations and data reduction 
are outlined in Sect. 3, while the data analysis is presented in 
Sect. 4. The morphology of the continuum stellar emission and 
ionized gas is presented in Sect. 5. Finally, the summary is given 
in Sect. 6. The second paper of the series will present the extinc- 
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Fig. 1. Luminosity distribution of the IRAS Revised BGS 
(Sanders et al. 2003, orange dashed histogram) compared with 
our sample of ULIRGs (black dashed histogram). The lu- 
minosities of the BGS galaxies have been re-calculated to 
consistently compare both datasets. The luminosity range is 
12.0<log(L /s /L o ) <12.6 with bins of 0.2. The two LIRGs of the 
sample have not been included in the graphic. 

tion of the sample. The third paper will focus on the ionization 
state, whereas the fourth will present the kinematics. Throughout 
the paper, we use £2a=0.7, £2^=0.3, and Ho=70 km s _1 Mpc -1 . 

2. Sample selection and properties 

The present sample of low-z ULIRGs was selected using the fol- 
lowing criteria: (1) to sample the IR luminosity range (Fig. [TJ; 
(2) to cover all types of nuclear activity, that is different exci- 
tation mechanisms such as Hn- (i.e., star formation activity), 
LINER- (i.e., superwinds, shocks), and Seyfert-like (i.e., pres- 
ence of an AGN); (3) to optimize the linear scales by selecting 
low-z galaxies; and (4) to span different phases of the interaction 
process. 

We selected our sample from those Sanders et al. (1988; 
1995), Melnick & Mirabel (1990), Leech et al. (1994), Kim et al. 
(1995), Lawrence et al. (1996), and Clements et al. (1996). Our 
final selection of 22 northern hemisphere ULIRGs (33 individ- 
ual galaxies) is shown in Fig. [2] and its main characteristics are 
presented in Table[T| In the following paragraphs, we describe 
the selection criteria in detail. 

The range of luminosities covered by the sample is 
1 1.8<log(L/#/L Q ) < 12.6, including two systems classified as 
LIRGs, but close enough to the low IR luminosity range of the 
ULIRG class. The luminosity distribution of the galaxies in our 
sample is similar to that of the complete flux-limitecQ IRAS 
Revised Bright Galaxy Sample (RBGS; Sanders et al.2003), 
which covers the entire sky surveyed by IRAS at Galactic lat- 
itudes |b|>5° (Fig. [T}. 

All the excitation mechanisms are well represented (see 
TablefT): 27%, 32%, and 18% of the ULIRGs are classi- 
fied as Hn, LINER, and Seyfert, respectively. The remain- 
ing 23% represent objects with unclear or mixed classifications. 
Using the 1 Jy sample, a flux-limited sample of ULIRGs iden- 
tified from the IRAS Faint Source Catalog, Kim et al. (1998) 

'threshold criteria S60/j m > 5.4 Jy, referred to the IRAS band centered 
in 60yum 
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derived the following distribution for ULIRGs: 28% for Hn, 
38% for LINER, and 34% for Seyfert. The H n and LINER frac- 
tions compare well with ours, whereas the Seyfert fraction is 
almost double. This difference is mainly caused by the wider lu- 
minosity covered by Kim et al. (1998; their sample goes up to 
L/fl=12.8L Q , whereas ours stops at 12.6L ), which results in a 
higher percentage of AGN (Veilleux et al. 1999). 

The galaxies are located at distances of between about 40 and 
900 Mpc, with an angular sampling that ranges from 0.2 to 3.1 
kpc arcsec -1 (see individual values in TablefTf. The median red- 
shift is 0.11, which compares well with 0.14, the median value 
of the 1 Jy sample (Kim & Sanders 1998). 

Finally, the ULIRGs in our sample exhibit a variety of mor- 
phologies, from interacting disks with a projected nuclear dis- 
tance up to 36 kpc, to close nuclei or individual nucleus sur- 
rounded by a common stellar envelope. There are detailed clas- 
sification schemes based on the interaction stage and morpho- 
logical features that lead to the existence of many classes and 
subclasses (see e.g., Surace 1998; Veilleux et al. 2002). We, how- 
ever, chose to separate our sample of ULIRGs into two broad 
categories using the projected nuclear distance as a discrimi- 
nator: (1) Pre-coalescence galaxies: those with a projected nu- 
clear separation larger than 1.5 kpc. (2) Post-coalescence galax- 
ies: those with a projected nuclear separation smaller or similar 
than 1 .5 kpc. This accurate distance calculation was possible be- 
cause of the high spatial resolution of the HST images (typical 
scale 07 1 per pixel). The selected distance evenly distributes the 
galaxies into the two categories, and is appropriate for separat- 
ing ULIRGs into early and late merger phases. In addition to 
this, theoretical models suggest that by the time the two nuclei 
have reached a separation of < 1 kpc, the stellar system has ba- 
sically achieved equilibrium, although their nuclei can still be 
separate structures (e.g., Mihos 1999; Mihos & Hernquist 1996; 
Bendo & Barnes 2000; Naab et al. 2006). However, it is impor- 
tant to clearly state that the use of projected nuclear distances 
may lead to some misclassification because galaxies identified 
as post-coalescence galaxies could well be in a pre-coalescence 
state. 

Considering all these aspects, we conclude that our sample 
is essentially representative of the ULIRG class. 



3. Observations and data reduction 

IFS observations were obtained with INTEGRAL, a fiber-based 
optical integral field system (Arribas et al. 1998) connected to 
the Wide Field Fibre Optic Spectrograph (WYFFOS; Bingham 
et al. 1994) and mounted on the 4.2 m William Herschel 
Telescope. The observations were carried out during a num- 
ber of observing runs between 1998 and 2004. We used three 
INTEGRAL configurations, namely the so-called standard bun- 
dles 1 (SB1), 2 (SB2), and 3 (SB3). The fiber diameters are 0745, 
079, and 2'.'7 for SB1, SB2, and SB3, respectively, resulting in 
the field of view (FoV) given in Table|2] The three bundles have 
a similar configuration in the focal plane. The majority of the 
fibers form a rectangular area centered on the object, whereas 
a subset of fibers that form an outer ring of 45" in radius is 
simultaneously used for measuring the sky (see configuration 
details in Table |2j. A remarkable capability of the instrument 
INTEGRAL is its flexibility in the bundle change, which allows 
us to select the most convenient instrument configuration de- 
pending on the seeing conditions. For the majority (~75%) of 
our cases, the SB2 bundle was the preferred one: the fiber size 



(079) is similar to the typical seeing of La Palma, and in general 
an entire ULIRG fits in the FoV (1670x1273). 

The spectra were taken with a 600 lines mirT 1 grating, which 
provides an effective spectral resolution (FWHM) of 6.0, 6.0, 
and 9.8 A for the SB1, SB2, and SB3 bundles, respectively 
The covered spectral range of interest was /1/14500-7000 A (rest- 
frame). The total integration time, fiber bundle used on each 
galaxy, and individual comments on some observations are pre- 
sented in Table |3] 

The data reduction was performed within the IRAF^] envi- 
ronment, followed the standard procedures applied to this type 
of data (see Arribas et al. 1997 and references therein), and can 
be summarized as follows. After subtracting an averaged BIAS 
frame, we proceeded to define the fibers apertures. This step in- 
volves the identification of information coming from each fiber, 
to preserve the full spatial and spectral information provided by 
IFS. For this reduction step, we used a sky flat image, which 
clearly indicates the position of spectra from the fibers along the 
detector because of its high signal over the entire wavelength 
range. The location of the INTEGRAL+ WYFFOS system at the 
Nasmyth focus of the telescope ensures its stability throughout 
a given night (del Burgo 2000), thus the information obtained 
for the sky flat can be used in all the scientific images observed 
during the same night. Once the fibers have been identified, we 
proceed to subtract the stray light and cross-talk contamination. 
Both signals create a background that has a diluting effect on the 
spectra, affecting in particular spectra from low surface bright- 
ness regions. The spurious light outside the defined apertures 
was analyzed, modeled two-dimensionally, and extracted from 
every image. The light transmitted by each fiber onto the detec- 
tor was afterwards measured, added along the spatial direction, 
and subtracted to obtain a final one dimensional spectrum for 
every single fiber. 

The one dimensional spectra were individually wavelength 
calibrated using a well-characterized lamp arc image. We 
checked the wavelength calibration using well identified sky 
lines, and obtained standard deviations of between 0.08 and 0. 18 
A. 

The next step in the data reduction involves the flat-field cor- 
rection. By combining the information provided by lamp and sky 
flats, which accounted for pixel-to-pixel and fiber-to-fiber varia- 
tions, respectively, we derived a response image, which was ap- 
plied to all science data. Finally, the sky contribution to the spec- 
tra was defined and subtracted using the combined information 
of the outer ring of fibers present in every bundle. We carefully 
checked that the fibers covering the sky were not contaminated 
by any contribution of the galaxy. 

For the flux calibration, several spectrophotometric standard 
stars^Jwere observed (BD+28 4211, Feige 34, BD +3326 42, 
and GD 153) using the same instrumental configuration and data 
reduction procedures as for the galaxies. Once the data reduc- 
tion was completed, we combined all the individual exposures 
(a minimum of three), improving the S/N and rejecting cosmic 



These values correspond to the old camera mounted on WYFFOS, 
which was used for the present observations. In August 2004, a new 
camera was commissioned for the instrument. See more details in 
http://www.iac.es/proyecto/integral 

^IRAF software is distributed by the National Optical Astronomy 
Observatory (NOAO), which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., in cooperation 
with the National Science Foundation. 

4 These stars are part of the Hubble Space Telescope spectrophoto- 
metric standards. 
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Fig. 2. HST (WFPC2/F814W) images of our sample of galaxies. From left to right, up to down, the projected nuclear separation 
decreases (i.e., the merger evolves from pre- to post-coalescence galaxies). Orientation is north up east to the left. The scale indicates 
5 arcsec. 
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.2. cont. Post-coalescence ULIRGs. 
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Table 1. Main characteristics of the sample of ULIRGs. The galaxies are sorted by decreasing nuclear distance. 



IRAS Name 



R.A.(I2000) 
(hh:mm:ss) 



DEC(J2000) 
(deg:mm:ss) 



Spectral Class 
(Optical) 



(Mpc) 



log(Lmf] 

(L G ) 



Scale 
(kpc arcsec~ 



Morpholog}^] 



IRAS 13156+0435 
IRAS 13342+3932 
IRAS 18580+6527 
IRAS 16007+3743 
IRAS 06268+3509 
IRAS 08572+3915 
IRAS 14348-1447 
Arp 299^ 

IRAS 13469+5833 
IRAS 12112+0305 
Mrk 463 

IRAS 06487+2208 
IRAS 11087+5351 
Mrk 273 
Arp 220 

IRAS 17208-0014 
IRAS 15250+3609 
IRAS 12490-1009 
IRAS 14060+2919 
IRASF 09427+1929 
IRAS 15206+3342 
Mrk 231 



13:18:07.6 
13:36:23.3 
18:58:06.3 
16:02:36.9 
06:30:12.6 
09:00:25.4 
14:37:38.4 
11:28:30.4 
13:48:40.3 
12:13:42.9 
13:56:02.8 
06:51:45.7 
11:11:36.4 
13:44:42.0 
15:34:57.1 
17:23:21.9 
15:26:59.4 
12:51:40.7 
14:08:18.9 
09:45:32.4 
15:22:38.0 
12:56:14.2 



+04:18:59.4 
+39:16:41.9 
+65:31:32.2 
+37:34:39.7 
+35:07:51.9 
+39:03:54.4 
-15:00:22.8 
+58:34:10.0 
+58:18:50.0 
+02:48:29.0 
+ 18:22:17.2 
+22:04:27.0 
+53:35:02.0 
+55:53:12.1 
+23:30:11.5 
-00:17:00.9 
+35:58:37.5 
-10:25:26.1 



+29 
+19 
+33 
+56 



04:46.9 
15:34.9 
31:35.9 
52:25.2 



H n/LINER (G) 
Seyl (E) 
Hn/Sey2 (G) 
LINER (H) 
Hn (Cfl 
LINER (Ef\ 
LINER (E) 
Hn/Sey2 (Dfl 
Hn (ifl 
LINER (E) 
Seyl/Sey2 (F, of] 
Hn(C) 
Seyl (G) 
LINER/Sey2 (B, 
LINER (E) 
LINER (AT 
LINER (HJ 
Hn/LINER( 
Hn(E) 
Hn (H) 
Hn(E) 
Seyl (E) 



0.113 
0.179 
0.176 
0.185 
0.169 
0.058 
0.083 
0.010 
0.158 
0.073 
0.050 
0.144 
0.143 
0.038 
0.018 
0.043 
0.055 
0.101 
0.117 
0.149 
0.124 
0.042 



525 

866 

850 

898 

813 

259 

378 

43 

755 

330 

222 

682 

677 

167 

78 

190 

245 

465 

545 

708 

580 

186 



12.13 
12.51 
12.26 
12.11 
12.5 
12.1 
12.39 
11.81 
12.31 
12.37 
11.81 
12.5 
12 
12.18 
12.20 
12.43 
12.09 
12.07 
12.18 
12.10 
12.27 
12.57 



il 



2.058 
3.026 
2.986 
3.100 
2.895 
1.130 
1.556 
0.205 
2.726 
1.395 
0.984 
2.522 
2.507 
0.749 
0.368 
0.844 
1.072 
1.854 
2.113 
2.600 
2.232 
0.832 



IP at 36.0 kpc 
IP at 25.1 kpc 
IP at 15.0 kpc 
IP at 14.2 kpc 
IP at 9.1 kpc 
IP at 6.1 kpc 
IP at 5.5 kpc 
IP at 5.0 kpc 
IP at 4.5 kpc 
IP at 4.0 kpc 
DN at 3.8 kpc 
DN at 1.5 kpc 
DNat 1.5 kpc 
DN at 0.7 kpc 
DN at 0.4 kpc 

SN 

SN 

SN 

SN 

SN 

SN 

SN 



"For the luminosity distances, the Wright (2006) cosmology calculator has been used. 

^(8- 1000 fim) was derived following Sanders & Mirabel 1996 (L(8 - lOOO^rn) = 4^D L 2 F«[L G ], with F F1R =\26 x 10~ 14 • (2.58/ 60 + 
fioo)[Wtn~ 2 ]. The quantities / 12 , /2s, feo, and f tm are the IRAS flux densities in ly at 12, 25, 60 and 100 /jm.). The IRAS fluxes were obtained from 
the IRAS Faint Galaxy Sample (Moshir et al. 1993, Vizier on line catalogue II/156A). 

TP means interacting pair, DN double nucleus and SN single nucleus. The distances in kpc represent the projected nuclear distance as measured 
in the F814WHST image. 

rf (C) Darling & Giovanelli 2006 

TRAS fluxes obtained from the IRAS point Source Catalogue. 
/(E) Kim, Veilleux & Sanders 1998. 

s Arp 299 is the system formed by NGC 3690 (to the east) and IC 694 (to the west). 
''(D) Garcfa-Marm et al. 2006. 
; Veilleux et al. 1999. 

'(F) Miller & Goodrich 1990; (G) Garcfa-Marin 2007. 
*IRAS fluxes obtained from the IRAS point Source Catalogue. 
'(B) Colina, Arribas & Borne 1999. 
'"(A) Arribas & Colina 2003. 
"(H) Veilleux et al. 1995. 



rays. Finally, when needed, a correction factor accounting for 
the percentage of flux not integrated by the central fiber (i.e., 
the fiber where the calibration star is centered) was used. Using 
the four spectrophotometric standard stars during six different 
observing runs, the average uncertainty obtained in the flux cal- 
ibration of the galaxies was in the range 10-15%. For further 
details about the flux calibration technique, we refer to Monreal- 
Ibero et al. (2007). 

We also retrieved broad-band imaging data for our galaxies 
from the HST archive. The images were taken with the Wide 
Field and Planetary Camera (WFPC2) using the F814W filter, 
and are available for the entire sample. This filter is equivalent to 
the ground-based Johnson-Cousins I (Origlia & Leitherer 2000). 
The images were calibrated on the fly, with the highest quality 
available reference files at the time of retrieval] 



5 More information about data process can be found at 
www.stsci.edu/hst/wfpc2/analysis/analysis.html 



4. Data analysis 

After data reduction, every galaxy has a set of spectra, each of 
them associated with a particular galaxy region observed by an 
individual fiber. The nuclear spectra of the systems under study 
are shown in Fig. [3] Depending on the instrument set-up, the an- 
gular scale of the spectra is 0745, 079, or 2'.'7 for the SB1, SB2, 
and SB3 bundles, respectively. In this context, nuclear spectra 
are the ones corresponding to the fiber closest to the continuum 
peak. 

Almost all galaxies exhibit nuclear emission, the most im- 
portant optical emission lines including those that are both 
strong (H/3, [Om]/U4959, 5007, Har+[Nn]Xl6548, 6484, and 
[Sn],M6716, 6731) and weak ([Oi]^6300). Some galaxies also 
exhibit other weak emission lines ([Ni]/15199, He 1/15876, and 
He 1/16678), as well as interstellar absorption lines (Nai/U5890, 
5896). 

To obtain maps of the most relevant optical lines tracing 
the ionization state, we first fitted the lines to Gaussian func- 
tions using the DIPSO package (Howarth & Murray 1988), in- 
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Fig. 3. Nuclear spectra of the sample of galaxies. The number of the INTEGRAL standard bundle used is specified next to the 
galaxy name. Each spectrum corresponds to a region of 0745, 079 and 2'.'7, depending on the bundle used (SB1, SB2, and SB3 
respectively). For Arp 220, we show the region of interest observed with the SB2 INTEGRAL bundle. For IRAS 13343+3932, we 
only show the western nucleus spectrum, since the eastern one was of low S/N ratio. The integration time for IRAS 13469+5833 
was insufficient to achieve a good S/N ratio. Its Ha+[Nn] complex was also strongly contaminated by sky lines. Hence, here we 
only show one spectrum as example. 



side the STARLINK environment It is well known that, under 
certain conditions, such as line blending, the fitting algorithms 
may derive different solutions. We therefore decided to apply 
some restrictions to constrain our results more reliably. The H/3 
and [Oi]/I6300 lines were individually fitted with a single com- 
ponent, with no particular constraint. We note that no Hy6 in ab- 
sorption was detected, and hence no correction was applied. The 
[Om]/L14959, 5007, and Ha+[Nn] complexes were fitted with 
two and three Gaussians, respectively, which were assumed to 
have identical kinematics. We applied these conditions based on 
the assumption that these line complexes trace identically the 
ionized gas kinematics. Additionally, we fixed the line intensity 
ratios of the [O in] and [N n] lines according to atomic parame- 
ters. For the [S n]AA6716, 6731 lines, we assume only that they 
share the same kinematics, since their line ratio variations can be 
used to trace the electron density. 

The spectral maps shown in this paper have all been derived 
using the values obtained with a single Gaussian component fit- 
ting for each line. In few cases, we detected additional kinemat- 
ical components in the line profiles, which may be indicative of 



the presence of winds, superwinds, and AGN flows. Some re- 
gions of a few galaxies are also consistent with the presence of 
a broad component in the hydrogen recombination lines. These 
particular cases will not be addressed here but in dedicated pa- 
pers about the kinematics and excitation conditions, because the 
purpose of the present work is to provide a general overview of 
the data. 

The mappings tracing the stellar component were derived us- 
ing a continuum emission line-free filter with a rectangular band- 
width of about 150 A for each galaxy (filters being centered ap- 
proximately at 4700 A and 6150 A restframe). 

Finally, the derivation of all these maps involves the ordering 
of the fibers according to their astronomical position (RA, Dec). 
To that aim, we used the free software IDA (Garcfa-Lorenzo et 
al. 2002), a tool specifically designed for the INTEGRAL IFU. 
To preserve the instrumental spatial resolution, the pixel size 
was equivalent to the fiber size used during the observations 
(see Figs. A.l - A. 28 i. The HST images and the spectral maps 



See http://www.starlink.rl.ac.uk/ 



were aligned by comparing with the peaks and external envelope 
structure of the red continuum map. 
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Table 3. Integral field spectroscopy observations log. The galaxies are sorted by decreasing nuclear distance. 



IRAS Name 


INTEGRAL 
(Bundle) 


INTEGRAL 
(FoV arcsec 2 ) 


Date 
month/year 


Exp. time 

(s) 


Comment on IFU data 


IRAS 13156+0435 


SB2, SB3 


32.9x25.3, 69.1x60.5 


03/02 


1800fF] 


[O in] non detected in N comp. 


IRAS 13342+3932 


SB3 


101.7x89.0 


03/02 


6xl50U 




IRAS 18580+6527 


SB2 


47.8x36.7 


04/01 


7x1500 




IRAS 16007+3743 


SB2 


49.6x38.1 


04/01 


6x1500 




IRAS 06268+3509 


SB2 


46.3x35.6 


03/02 


6x1500 




JTI 

IRAS 08572+3913^ 


SB2 


18.0x13.9 


04/99 


9x1800 




IRAS 14348-1447^ 


SB2 


1.3x1.0 


04/98 


4x1800 




1 — 1 

Arp 299H 


SB2, SB3 


3.3x2.5, 6.9x6.0 


01/04 


1440f|] 




IRAS 13469+5833 


SB1 


21.3x17.4 


03/02 


5xl80U 


Sky line to Ha+fNii] complex 












Not enough S/N 


IRAS 12112+0305 fl F] 


SB2 


27.9x19.5 


04/98 


5x1800, 4x1500 




Mrk 463 


SB2 


15.7x12.1 


04/01 


5x900 




IRAS 06487+2208 


SB1 


19.7x16.1 


03/02 


5x1500 




IRAS 11087+5351 


SB1 


19.5x16.0 


03/02 


6x1500 




Mrk 273" 


SB2 


12.0x9.2 


04/98 


3x1500 




Arp 220" 


SB2, SB3 


5.9x4.5, 12.4x10.8 


05/00 


4800(Q 


[O m] not detected, FI/3 not covered 


IRAS 17208-0014° 


SB2 


13.5x10.4 


04/99 


4x1950 




IRAS 15250+3609" 


SB2 


17.1x13.2 


04/98 


5x1800 




IRAS 12490-1009 


SB2 


29.7x22.8 


03/02 


4x1500 




IRAS 14060+2919 


SB2 


1.9x1.4 


03/02 


4x1800 




IRAS 09427+1929 


SB3 


87.4x76.4 


03/02 


232(F] 


Short inhomogeneous exposures 


IRAS 15206+3342° 


SB2 


35.7x27.4 


04/99 


4x1800 




Mrk 231 


SB1 


6.5x5.3 


03/02 


4x900 


H/? not covered 



Broad emission lines 



"3x1500 s for the South component with the SB2 bundle; 3x1800 for the North component with the SB2 bundle, 5x1500 s for the entire 
system with the SB3 bundle. 

''Previously published galaxy. See Colina, Arribas & Monreal-Ibero 2005, and references therein. 
' Garcfa-Marin et al. 2006. 

rf 3xl200 s for each individual galaxy of the system, NGC 3690 and IC 694 observed with the SB2 bundle. 3x1800 s for the entire system as 
observed with the SB3 bundle. 

''For IRAS 121 12+0305 we present the results obtained by combining two different pointings of the SB2 INTEGRAL bundle. 
^6x1500 s for the nuclear region observed with the SB2 bundle, 26x1500 s for the nucleus and the extended nebula as observed with the SB3 
bundle. 

*1500 s, 520 s, 300 s. 



5. Observed morphologies of the stellar and 
ionized gas components 

5. 1 . General trends 

We present and discuss the differences between the structure 
of the stellar and ionized gas components for our sample of 
ULIRGs. 

The high resolution HST images provide evidence for very 
complex stellar structures, with bright nuclear regions, dust 
structures, and knots of star formation. In the pre-coalescence 
systems, the large scale structure is characterized by merger- 
induced structures, i.e., bridges, tails, and plumes connect- 
ing the individual galaxies. The original spiral arms are dis- 
rupted and form kpc-size tidal tails with, in some cases, chains 
of star-forming regions (Monreal-Ibero et al. 2007). For the 
post-coalescence galaxies, the large scale structure traces an 
outer common envelope, which resembles an elliptical structure 
(Garcfa-Marin 2007), in addition to late merger features such as 
disordered nuclear regions and tails. 

In spite of the different angular sampling, the INTEGRAL- 
based continuum images of the galaxies obtained at wavelengths 
blueward of Ha and Hfi recover the stellar structure observed 
in the high angular resolution HST images. This occurs even in 
the low-surface brightness regions, reassuring us that all the ob- 



served stellar and ionized gas structures are real. Only small dif- 
ferences between the INTEGRAL red and blue continua are de- 
tected (e.g., Arp 299, IRAS 08572+3915, IRAS 12112+0305), 
and they can probably be explained by extinction effects. 

The general morphology of the ionized gas low surface 
brightness regions is similar to that of the stellar component, al- 
though in about 25% of the galaxies structural variations such 
us differences in shape or the presence of external clumps are 
associated with gas in high excitation states (Arp 299, Mrk 463, 
Mrk 273, IRAS 18580+6527, and IRAS 11087+5351). These 
differences are mainly traced by the [Oni]/15007. In the pre- 
coalescence systems, one of the optical nuclei is usually a weak 
line emitter, with a surface brightness similar to that of the inter- 
face between the galaxies. 

The largest structural variations between the stellar and ion- 
ized gas components are observed in the high surface brightness 
regions (i.e., the main body of the galaxies), and are detected 
as differences in the location of the emission peaks (see the dis- 
tance between continuum and emission lines peaks for individ- 
ual ULIRGs in Table|4]l. In what follows, we consider a signifi- 
cant shift as those equal to or larger than 1 kpc (within the errors 
given in TableQ. 

There is also evidence of substructures in the different ion- 
ization states of the gas. These variations may be due to mecha- 
nisms such as star formation, AGNs, and shocks caused by tidal 
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effects. To trace this more reliably, we have selected respectively 
the strong hydrogen recombination line Ha, the [Om]/I5007 
line, which traces hard ionizing photons, and [Oi]/I6300, which 
is considered to be a sensitive shock tracer. Other effects, such as 
extinction and differences in stellar populations and in metallici- 
ties, may also play a role, but the detailed analysis of the excita- 
tion mechanisms in ULIRGs is beyond the scope of the present 
paper. In the following, we present the main structural character- 
istics of the pre- and post-coalescence systems. 

Some ULIRGs of the sample were not included for a variety 
of reasons, such as the lack of [Om]/I5007 detection (Arp 220), 
the low S/N of their data precluding the detection of emission 
lines (e.g., IRAS 13469+5833), poor spatial resolution (e.g., 
IRAS 13342+3932), and severe AGN contamination of the host 
galaxy light that makes it impossible to distinguish the ionized 
gas component (Mrk231, see Appendix). The analysis presented 
here was carried out for 18 ULIRGs (i.e., 20 INTEGRAL point- 
ings, see Table[3]). This includes two individual pointings for the 
spatially separated galaxies that form the systems Arp 299 and 
IRAS 13156+0435. On average, the linear scale corresponding 
to our angular resolution is about 1.8 kpc. Therefore, all the re- 
sults presented here should be considered valid on this scale. 

5.2. Pre-coalescence ULIRGs 

The ionized gas traced by the Ha maps in the pre-coalescence 
systems exhibit emission on scales of about 5-10 kpc measured 
from the stellar nuclei (see Figs, A.l - A. 16 1. Their intensity 



peak are located at different positions: in the brightest stellar nu- 
clei (e.g., IRAS 18580+6527), in the main body/overlapping re- 
gion between the galaxies (e.g., NGC 3690), in one of the tails 
(e.g., IRAS 16007+3743), in an extra-nuclear stellar knot (e.g., 
IRAS 121 12+0305), or in the secondary nucleus (e.g. Mrk 463). 
In about 65% of the galaxies studied, there are significant off- 
sets between the Ha and the stellar continuum peaks, with 
projected separations of between 1-8 kpc. Extra-nuclear star- 
forming regions, mainly located in tidal tails, are also detected 
as secondary Ha peaks (IRAS 16007+3743, IRAS 08572+3915, 
IRAS 14348-1447, and IRAS 12112+0305). These additional 
light peaks are identified as Tidal Dwarf Galaxies (TDGs), can- 
didates (Monreal-Ibero et al. 2007). The H/3 structure is rather 
similar to that of Ha, but it is intrinsically weaker and more af- 
fected by extinction, especially in the galaxy nuclei. 

The overall morphology, relative brightness of the nu- 
clei, and the peak distribution of the brightest forbidden lines 
[Nn]/16484, and [S n]/U6716, 6731 are in good agreement with 
that of Ha. Nonetheless, there are some cases where small scale 
differences in their peak positions are found (IRAS 18580+6527 
and IRAS 12112+0305), and one galaxy where their peaks are 
located 8.4 kpc away from the one of Ha (IRAS 16007+3743). 
Given the wavelength proximity of these emission lines, it is un- 
likely that these differences are produced by any extinction ef- 
fects. They are probably produced instead by the presence of 
different ionization mechanisms associated with tidal shocks. 

The results for the highly excited line [Om]/I5007 are an 
interesting case study. As for Ha, in about 50% of the pre- 
coalescence systems, there are significant shifts between the stel- 
lar continuum and the nebular emission peaks. In comparison 
with other emission-line maps, we find that the [O m]/l5007 peak 
coincides with that of Ha in 70% of the studied cases. In some 
galaxies, there are regions that appear to be dominated by lo- 
cal ionization sources. For instance, in IRAS 08572+3915 and 
IRAS 12112+0305, the [Om]^5007 line has its maxima in re- 
gions identified as TDG candidates in the stellar continuum im- 



ages. These TDG candidates are also secondary Ha emitters. 
These differences may be caused by the effects of extinction 
on [Om]/15007, which are more important in the nuclei than in 
these external regions. In the case of IC 694, the peak of the 
[Om]/I5007 coincides with a region of the galaxy believed to 
be one of the original spiral arms disrupted by the merger. It is 
therefore also bright in the stellar component and in Ha. 

The galaxy IRAS 18580+6527 also deserves special men- 
tion, since the extra-nuclear local maxima located at about 



10 kpc from the eastern nucleus (also bright in Ha, see Fig. A. 5 
Aa — 4, A6~-3) is not associated with any particular stellar 
mass concentration, but is rather dominated by nebular emis- 
sion. Interestingly enough, this galaxy is classified as a Seyfert. 
As explained below, similar extra-nuclear highly excited nebu- 
lae are detected in the post-coalescence ULIRGs Mrk 273 and 
IRAS 11087+5351. 

Approximately 40% of the pre-coalescence ULIRGs present 
shifts of up to 6 kpc between the peaks of the red contin- 
uum and the [Oi]/I6300 line. However, in some cases such as 
IRAS 06268+3509 and Mrk 463, the line maxima is associated 
with the stellar secondary nucleus of the system. This would 
mean that in general the [Oi]/I6300 activity is mainly related 
to the nuclei. 



5.3. Post-coalescence ULIRGs 

For the morphologically evolved ULIRGs, the comparison be- 
tween the stellar and ionized gas components provides different 
results from above. The most important one is that all the emis- 
sion tends to be concentrated in the nuclear region. The pre- 
viously reported structural differences between the stellar and 
ionized gas are not so common for this type of ULIRGs. The 
spectral maps indicate that in all cases (see Figs. |A. 17] - |A.28| > 
the Ha peak coincides with the peak of the stellar emission (i.e. 
within the central kpc), indicating a nuclear concentration of the 
excitation sources. The Ha overall structure generally coincides 
with that of the red continuum. One galaxy (IRAS 15250+3609) 
contains a TDG candidate that is also an Ha emitter (Monreal- 
Ibero 2007). For Mrk 273, Arp 220, and IRAS 11087+5351, 
there are also secondary extra-nuclear structures that do not cor- 
relate spatially with the stellar distribution. Likewise for the pre- 
coalescence systems, the Ha map generally agree with those of 
H/3, although the latter is more affected by extinction that can 
cause their maxima position not to be coincident. 

The behavior of the [Nn]^6484, and [S n]/U6716, 6731 fol- 
lows that of Ha, and only small scale variations are measured in 
IRAS 15250+3609. 

As in the pre-coalescence systems, there are two galaxies 
(IRAS 11087+5351 and Mrk 273) that exhibit a peak in the 
[Om]/I5007 emission within the extra-nuclear regions (7.2 and 
5.1 kpc, respectively). These regions are not dominated by lo- 
cal ionization sources, because there is no particular stellar mass 
concentration there. Interestingly, both galaxies have a Seyfert 
classification (see Table[T]l. Secondary Ha structures are also de- 
tected in these extra-nuclear regions. 

The [Oi]/I6300 emission line maps also follow the nuclear 
concentration trend exhibited by the post-coalescence systems, 
with only the maxima position of IRAS 11087+5351 being 
shifted with respect to that of the continuum and almost coin- 
ciding with the [Om]/I5007 peak. 
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6. Summary 

This is the first paper in a series to study in detail the two- 
dimensional morphology, dust distribution, excitation processes, 
and kinematics of a representative sample of 22 local ULIRGs. 
To achieve this, we have used optical IFS data combined with 
high resolution HST images. In this paper, we have described the 
sample selection and the observations, and presented the mor- 
phological properties of the continuum and ionized gas compo- 
nents. The ULIRGs were selected to cover a wider range of the 
IR luminosity distribution, different merger phases, and to ex- 
hibit a variety of nuclear activities (Hn, LINER, Seyfert). Using 
the projected nuclear distance as a criteria, we have adopted a 
simple classification scheme that divides the sample into pre- 
and post-coalescence galaxies, and identifies the different be- 
haviors at different stages of the merger. The main results of this 
paper can be summarized as: 

- Despite the different resolution, the structure of the 
INTEGRAL red continua is consistent with that of the HST 
F814W band. The red and blue continua may exhibit slight 
differences caused by the effect of the extinction. In contrast, 
the structure of the warm ionized gas can be significantly de- 
coupled from that of the stellar continuum. These variations 
are due to the different spatial distribution of the ionization 
sources present in the galaxy structure, and to the dust distri- 
bution. 

- We have compared the stellar and ionized gas structures 
of ULIRGs undergoing different merger phases (pre- and 
post-coalescence), and found structural variations between 
the gaseous and the stellar components, with offsets in the 
emission peak positions of up to about 8 kpc in the pre- 
coalescence systems. For 64% of ULIRGs in an early inter- 
action phase, the Ha peak does not coincide with the stellar 
maxima. In contrast, galaxies undergoing advanced mergers 
have their Ha peak located in the nuclear regions. 

- Four pre- and one post-coalescence ULIRGs appear to con- 
tain TDGs candidates, which we identify with Ha emission 
associated with stellar emission knots. 

- The ionization structure traced by different emission lines 
has been also studied. We detected variations in lines with 
similar rest-frame, meaning that the differences are unlikely 
to be explained by the presence of dust. This is more com- 
mon in the pre-coalescence ULIRGs, and it is probably a 
consequence of different excitation mechanisms. Differences 
in metallicity and stellar populations may also be playing a 
role. 

- The [Om] line traces highly excited extra-nuclear clouds 
in both pre- and post-coalescence systems, with no rele- 
vant stellar counterpart. These regions are also secondary Ha 
emitters. All these galaxies have a Seyfert nuclei, classified 
using optical emission lines ratios. 

- The peak of the [O i] line is shifted with respect to the stellar 
one, but in general its emission is nuclear. 

- This analysis infers that there is an evolutionary trend in the 
ionized gas behavior. In the pre-coalescence ULIRGs, the 
structural differences are remarkable, and the extra-nuclear 
regions, out to distances of about 8 kpc, play an important 
role and contribute to the excitation level. In contrast, the 
post-coalescence ULIRGs tend to have their activity concen- 
trated in the nucleus. The only exception to this are the extra- 
nuclear nebulae traced by the [Om]. We have detected them 
in both galaxy types, and they appear to be related to the 
presence of a Seyfert nucleus. 



Acknowledgements. We would like to thank Almudena Alonso-Herrero 
for useful comments and discussions. This paper uses the plotting 
package jmaplot, developed by Jesus Maiz-Apellaniz. http://dae45. 
iaa.csic.es:8()80/~jmaiz/software. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 

This work has been supported by the Spanish Ministry of Education and 
Science, under grant BES-2003-0852, project AYA2002-01055. MG-M is sup- 
ported by the German federal department for education and research (BMBF) 
under the project numbers: 50OS0502 & 50OS0801. AMI is supported by 
the Spanish Ministry of Science and Innovation (MICINN) under the program 
"Specialization in International Organisms" ref. ES2006-0003. 



References 



Alonso-Herrero, A. 2009, A&A submitted 

Arribas, S., Colina, L., Monreal-Ibero, A., Alfonso, J., Garcia-Marin, M. & 
Alonso-Herrero, A. 2008, A&A, 479, 687 

Arribas, S., Colina, L., 2003, ASP Conference Proceedings, 297, 24-28 
Arribas, S., Colina, L. & Borne, K. 1999, Ap&SS, 266, 143 
Arribas, S. et al. 1998, Proc. SPIE, 3355, 821 

Arribas, S., Mediavilla, E., Garcia-Lorenzo, B., & del Burgo, C. 1997, ApJ, 
490, 227 

Bendo, G. J. & Barnes, J, E. 2000, MNRAS, 316, 315 

Bingham, R. G., Gellatly, D. W., Jenkins, C. R., & Worswick, S. P. 1994, Proc. 
SPIE, 2198, 56 

Borne, K.D., Bushouse, H, Lucas, R. A., & Colina, L. 2000, ApJ, 529L, 77 
Bushouse, H. A., Borne, K. D., Colina, L., Lucas, R. A., Rowan-Robinson, M., 
Baker, A. C, Clements, D. L., Lawrence, A. & Oliver, S. 2002, ApJS, 138, 1 
Caputi, K. I. et al. 2007, ApJ, 660, 97 
Caputi, K. I. et al. 2006a ApJ, 637, 727 

Caputi, K. I. McLure, R. J., Dunlop, J. S., Cirasuolo, M. & Schael, A. M. 
2006b, MNRAS, 366, 609 

Clements, D. L., Sutherland, W. J., Saunders, W., Efstathiou, G. P., McMahon, 
R. G, Maddox, S., Lawrence, A. & Rowan-Robinson, M. 1996, MNRAS, 279, 
459 

Colina, L., Arribas, S. & Monreal, A. 2005, ApJ, 621, 725 
Colina, L., Arribas, S. & Borne, K. D. 1999, ApJ, 527L, 13 
Darling, J. & Giovanelli, R. 2006, AJ, 132, 2596 

Dasyra, K. M., Tacconi, L. J., Davies, R. I., Naab, T, Genzel, R., Lutz, D., 
Sturm, E., Baker, A. J., Veilleux, S., Sanders, D. B., & Burkert, A. 2006, ApJ, 
651, 835 

Elbaz, D., Cesarsky, C. J., Chanial, P., Aussel, H, Franceschini, A., Fadda, D. 
& Chary, R. R. 2002, A&A, 384, 848 

Farrah, D., Bernard-Salas, J., Spoon, H. W. W., Soifer, B. T., Armus, L., Brandl, 
B., Charmandaris, V., Desai, V., Higdon, S., Devost, D. & Houck, J. 2007, ApJ, 
667, 149 

Farrah, D., Rowan-Robinson, M., Oliver, S., Serjeant, S., Borne, K., Lawrence, 
A., Lucas, R. A., Bushouse, H, & Colina, L. 2001, MNRAS, 326, 1333 
Forster Schreiber, N. M. et al. 2009, arXiv0903.1872 
Forster Schreiber, N. M. et al. 2006, ApJ, 645, 1062 

Garcia-Lorenzo, B., Acosta-Pulido, J. A., & Megias-Fernandez, E. 2002, ASP 
Conference Proceedings, 282, 501 

Garcia-Marin , M., Colina, L., Arribas, S., Alonso-Herrero, A. & Mediavilla, E. 
2006, ApJ, 650, 850 

Garcia-Marin 2007, PhD, Universidad Autonoma de Madrid 
Genzel, R. et al. 2006, Nature, 442, 786 

Genzel, R., Tacconi, L. J., Rigopoulou, D., Lutz, D., & Tecza, M. 2001, ApJ, 
563,527 

Genzel, R., Lutz, D., Sturm, E., Egami, E., Kunze, D., Moorwood, A. F. M., 
Rigopoulou, D., Spoon, H. W. W., Sternberg, A., Tacconi-Garman, L. E., et 
al. 1998, ApJ, 498, 579 

Howarth, I. D., & Murray, J. 1988, DIPSO A Friendly Spectrum Analysis 
Program (Starlink User Note 50; Chilton: Rutherford Appleton Lab.) 
Imanishi, M., Dudley, C. C, Maiolino, R., Maloney, P. R., Nakagawa, T. & 
Risaliti, G. 2007, ApJS, 171,72 

Kim, D. C, Veilleux, S. & Sanders, D. B. 1998, ApJ, 508, 627 

Kim, D.-C, Sanders, D. B., Veilleux, S., Mazzarella, J. M., & Soifer, B. T. 

1995, ApJS, 98, 129 

Lagache, G. Puget, J-L. & Dole, H. 2005 ARA&A, 43, 727 
Le Floc'h, E. et al. 2005, ApJ, 632, 169 
Le Floc'h, E. et al. 2004, ApJS, 154, 170 

Leech, K. J., Rowan-Robinson, M., Lawrence, A. & Hughes, J. D. 1994, 
MNRAS, 267, 253 



Garcfa-Marin et al.: IFS of a Representative Sample of ULIRGs. I. The Data 



11 



Lonsdale, C, Farrah, D. & Smith, H. 2006, Astrophysics Update 2, edited 

by John W. Mason. ISBN 3-540-30312-X. Published by Springer Verlag, 

Heidelberg, Germany, 285 

Lonsdale, C. et al. 2004, ApJS, 154, 54 

Low, J., Kleinmann, & D. E. 1968, AJ, 73, 868 

Melnick, J. & Mirabel, I. F. 1990, A&A, 23 1L, 19 

Mihos, J. C. & Hernquist, L. 1996, ApJ, 464, 641 

Mihos, J. C. 1999, Ap&SS, 266, 195 

Miller, J. S. & Goodrich, R. W. 1990, ApJ, 355, 456 

Monreal-Ibero, A., Colina, L., Arribas, S. & Garci'a-Marin, M. 2007, A&A, 
472, 421 

Monreal-Ibero, A., Arribas, S. & Colina, L. 2006, ApJ, 637, 138 
Monreal-Ibero, A. 2004, PhD Thesis, Universidad de La Laguna 
Moshir, M., Copan, G., Conrow, T., McCallon, H., Hacking, P., Gregorich, D., 
Rohrbach, G., Melnyk, M., Rice, W., & Fullmer, L. 1993, VizieR On-line Data 
Catalog: II/156A 

Naab, T., Jesseit, R. & Burkert, A. 2006, MNRAS, 372, 839 

Nardini, E., Risaliti, G., Salvati, M., Sani, E., Imanishi, M., Marconi, A. & 

Maiolino, R. 2008, MNRAS, 385L, 130 

Origlia L., & Leitherer C. 2000, AJ, 1 19, 2018 

Perez-Gonzalez, P. G., Rieke, G. H., Egami, E., Alonso-Herrero, A., Dole, H., 
Papovich, C, Blaylock, M., Jones, J., Rieke, M., Rigby, J. & 4 coauthors 2005, 
ApJ, 630, 82 

Rieke, G. H. & Low, F. J. 1972, ApJ, 176L, 95 

Risaliti, G., Maiolino, R., Marconi, A., Sani, E., Berta, S., Braito, V., Ceca, R. 

D., Franceschini, A., & Salvati, M. 2006, MNRAS, 365, 303 

Rupke, D. S., Veilleux, S. & Sanders, D. B. 2005, ApJ, 632, 751 

Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, J. A. & Soifer, B. T. 

2003, AJ, 126, 1607 

Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749 

Sanders, D. B., Soifer, B. T., Elias, J. H., Neugebauer, G., & Matthews, K. 
1988, ApJ, 328L, 35 

Scoville, N. Z., Evans, A. S., Thompson, R., Rieke, M., Hines, D. C, Low, F. 

J., Dinshaw, N., Surace, J. A. & Armus, L. 2000, AJ, 1 19, 991 

Smail, Ian, Ivison, R. J., & Blain, A. W. 1997, ApJ, 490L, 5 

Schmitt, H. R., Donley, J. L., Antonucci, R. R. J., Hutchings, J. B., Kinney, A. 

L. & Pringle, J. E. 2003, ApJ, 597, 768 

Soifer, B. T, Rowan-Robinson, M., Houck, J. R., de Jong, T., Neugebauer, G., 
Aumann, H. H., Beichman, C. A., Boggess, N., Clegg, P. E., Emerson, J. P. & 
6 coauthors 1984, AjP, 278L, 71 

Surace, J. A., Sanders, D. B., Vacca, W. D., Veilleux, S. & Mazzarella, J. M. 
1998, ApJ, 492, 116 

Tacconi, L. J., Genzel, R., Lutz, D., Rigopoulou, D., Baker, A. J., Iserlohe, C, 
& Tecza, M. 2002, ApJ, 580, 73 

Veilleux, S., Kim, D.C., & Sanders, D. B. 2002, ApJS, 143, 315 

Veilleux, S., Kim, D.C., & Sanders, D. B. 1999, ApJ, 522, 113 

Veilleux, S., Kim, D.C., Sanders, D. B., Mazzarella, J. M., & Soifer, B. T. 1995, 

ApJS, 98, 171 

Yan, H. et al. 2004, ApJ, 616, 63 



12 Garcfa-Marin et al.: IFS of a Representative Sample of ULIRGs. I. The Data 

Table 4. Separation between the red continuum and the emission line gas peaks. 



Galaxy 


Nuc Sep 


[Oi] 


Ho- 


[Om] 


HyS 


[Sn] 


[Nil] 


Comment 




(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 





IRAS 13156+0435N 36.0 1.8±0.9 1.8±0.9 2.0±0.9 1.8±0.9 1.8±0.9* 1.8±0.9 

IRAS 13156+0435S 36.0 1.8±0.9 1.8±0.9 1.8+0.9 1.8±0.9 1.8±0.9 1.8±0.9 

IRAS 18580+6527 15.0 0.0+1.2 2.7+1.2 2.7+1.2 2.7+1.2 0.0+1.2 2.7+1.2 

IRAS 16007+3743 14.2 3.9+1.4 8.4+1.4 8.4+1.4 8.4+1.4 0.0+1.4 0.0+1.4 

IRAS 06268+3509 9.1 5.8+1.3 2.6+1.3 2.6+1.3 3.7+1.3 2.6+1.3 2.6+1.3 

IRAS 08572+3915 6.1 5.9+0.5 5.9+0.5 5.1+0.5 1.0+0.5 5.9+0.5 5.9+0.5 

IRAS 14348-1447 5.5 0.0+0.7 0.0+0.7 0.0+0.7 0.0+0.7 0.0+0.7 0.0+0.7 

Arp 299/NGC 3690 5.0 0.4+0.1 1.3+0.1 1.2+0.1 1.2+0.1 1.2+0.1 1.3+0.1 

Arp299/IC694 5.0 0.8+0.1 0.9+0.1 1.3+0.1 1.3+0.1 0.9+0.1 0.9+0.1 

IRAS 12112+0305 4.0 1.3+0.6 2.0+0.6 1.2+0.6 2.5+0.6 1.2+0.6 1.2+0.6 

Mrk463 3.8 3.5+0.4 3.5+0.4 3.5+0.4 3.5+0.4 3.5+0.4 3.5+0.4 

IRAS 06487+2208 1.5 1.1+0.6 1.1+0.6 1.1+0.6 0.0+0.6 0.0+0.6 1.1+0.6 

IRAS 11087+5351 1.5 6.6+0.6 1.6+0.6 7.2+0.6 7.2+0.6 1.6+0.6 1.6+0.6 

Mrk273 0.7 0.7+0.3 0.7+0.3 5.1+0.3 5.1+0.3 0.7+0.3 0.7+0.3 

Arp220 0.4 0.0+0.2 0.5+0.2 N/A N/A 0.3+0.2 0.3+0.2 

IRAS 12490-1009 0.0 1.7+0.8 1.7+0.8 0.0+0.8 0.0+0.8 1.7+0.8 1.7+0.8 

IRAS 14060+2919 0.0 1.9+0.9 1.9+0.9 1.9+0.9 1.9+0.9 1.9+0.9 1.9+0.9 

IRAS 15206+3342 0.0 0.0+1.0 0.0+1.0 0.0+1.0 0.0+1.0 0.0+1.0 0.0+1.0 

IRAS 15250+3609 0.0 0.0+0.5 1.4+0.5 1.4+0.5 1.4+0.5 0.0+0.5 0.0+0.5 

IRAS 17208-0014 0.0 0.0+0.4 0.0+0.4 1.1+0.4 1.1+0.4 0.0+0.4 0.0+0.4 



Activity concentrated in the nucleus 
Activity concentrated in the nucleus 
Secondary [O in] structure. 
Ha, [Oi] and [Om] offsets. TDG candidates. 
[O i] peaks on secondary nucleus. 
Ha peaks in secondary nucleus. 
Extended [O in] nebula. TDG candidates. 
TDG candidates. 
Ha, [Oi] and [O m] offset. 
Ha, [O i] and [O ra] offset. 
Ha and [Om] offsets. TDG candidates. 
Ionized gas peaks in secondary nucleus. 
Activity concentrated in the nucleus 
Extended [O m] nebula. 
Extended [O in] nebula. 

Activity concentrated in the nucleus 
Activity concentrated in the nucleus 
Activity concentrated in the nucleus 
Activity concentrated in the nucleus.TDG candidate. 
Activity concentrated in the nucleus. 
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Appendix A: Stellar component and ionized gas 
maps for the sample of ULIRGs: 

In this Appendix, the maps of interest for all the galaxies of the 
sample, as obtained with the INTEGRAL system, are presented. 
The emission-line free stellar continuum, along with the most 
relevant optical-emission line (Hfi, [Om]/I5007, [Oi]/i6300, 
Ha+[Nn]X16548, 6584, [S n],U)6716, 6731) maps are shown. 
Complementary optical (0.8 /im) HST images are presented too. 
The galaxies are sorted by decreasing nuclear separation. 
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Fig.A.l. Emission-line free stellar continua and emission-line images of IRAS3 156+0435 obtained with the INTEGRAL SB3 
bundle. Contours represent the red continuum. The HST-WFPC2 F814W image is shown for comparison. All the images are shown 
in a logarithmic scale. For a better comparison, the emission line maps have the same logarithmic color scale. The color code is 
given in relative flux units. North is up, east to the left. The scale represents 5 kpc. The arrows indicate the nuclei. 
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Fig. A.2. As Fig A. 1 but for a galaxy observed with the SB2 bundle. 
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Fig. A.3. As Fig A. 1 but for a galaxy observed with the SB2 bundle. 
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Fig. A.4. As Fig 



A.l 
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Fig. A.7. As Fig A. 1 but for a galaxy observed with the SB2 bundle. 
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Fig. A.8. As Fig A. 1 but for a galaxy observed with the SB2 bundle. 
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Fig. A.9. As Fig A. 1 but for a galaxy observed with the SB2 bundle. This pointing covers the southern tail of IRAS 
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Fig. A. 11. As FigA.l In the Arp 299 system, NGC 3690 is located to the west, and IC 694 to the east. In this case the scale 
represents 2 kpc. 
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Fig. A.12. As Fig A.l but for a galaxy observed with the SB2 bundle. In this case the scale represents 2 kpc. 
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Fig. A.14. As Fig A.l but for a galaxy observed with the SB1 bundle. In this case, the Hq-+[Nii] complex was strongly contaminated 
by a sky line. Also, the integration time was not sufficiente to achieve enough S/N ratio in the majority of the emission lines. 
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Fig. A.15. As Fig A.l but for a galaxy where we have combined two pointings made with the SB2 bundle. 
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Fig.A.21. As Fig 



A.l 



but for a mosaic obtained combining three pointings with the SB3 bundle. 
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Fig. A.22. As Fig A. 1 but for a galaxy observed with the SB2 bundle. 
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Fig. A.23. As Fig A. 1 but for a galaxy observed with the SB2 bundle. 
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Fig. A.26. As Fig A.l but for a galaxy observed with the SB3 bundle. 
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Fig. A.27. As Fig A.l but for a galaxy observed with the SB3 bundle. 
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Fig. A.28. As Fig A. 1 but for a galaxy observed with the SB 1 bundle. The strong contamination coming from the Seyfert nucleus 
precludes the lines detection. Only the Ha+[Nn] complex has being detected as a single structure. 
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Table 2. Characteristics of the INTEGRAL fiber bundles 



Bundle Fiber core (0 


a 


# Fiberf 


Fo\ 


C 


Sky-rin^ d 


SBl 0.45 ~ 205(175+7 
SB2 0.90 219(189+3 
SB3 2.70 135(115+2 


D) 7.80x6740 90 
3) 16.0x12.3 90 
D) 33.6x29.4 90 



"Fiber core diameter (arcsec). 
* Total number of fibers (rectangle+sky ring). 
c Spatial coverage of the central rectangle in arcsec. 
''Diameter of the external sky ring in arcsec. 



